Exact boson mapping of the reduced BCS pairing Hamiltonian 
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Abstract 

An exact boson mapping of the reduced BCS (equal strength) pairing Hamiltonian is considered. 
In the mapping, fermion pair operators are mapped exactly to the corresponding bosons. The image 
of the mapping results in a Bose-Hubbard model with level dependent hopping. Though the resul- 
tant Bose-Hubbard Hamiltonian is non-Hermitian, all eigenvalues are real when Uk/t < 1, where k 
is the total number of bosons. When U/t = 1, a part of spectrum of the Bose-Hubbard Hamiltonian 
corresponds exactly to the whole spectrum of the reduced BCS pairing Hamiltonian. 
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Pairing is one of important residue interactions in many areas of physics, especially in the stud- 
ies of superconductors, I 1 ' 2 ' nuclear systems, ^ metallic clusters, I 4, 5 ' and liquids.' 6 ' The limitations of 
£> ' the Bardeen-Cooper-Schrieffer (BCS) and Hartree-Fock-Bogolyubov (HFB) methods^ 1,7 ' for finding ap- 

proximate solutions of finite nuclear systems and nanoscale metallic grains are well understood. I 4, 8 1 
Fortunately, the reduced BCS (equal strength) pairing model was proved to be exactly solvable following 
Richardson's early work! 9, 10 1 and studies based on the Gaudin algebraic Bethe ansatz method, I 11 ' which 
f*^ ' has received a lot of attention recently. I 12, 13 ' 

. On the other hand, there is also a long history in search for appropriate boson mapping methods or 

boson expansions for nuclear many-body systems^ 14 ! which can also be carried out for other many-fermion 
systems. It is well known that in most circumstances pairs of fermions exhibit boson-like behavior. In 
such approaches, the degrees of freedom of fermion pairs are directly replaced by exact boson degrees of 
freedom. These methods are potentially helpful in describing collective motion in terms of boson degrees 
of freedom to avoid usual difficult fermionic formulation since boson operators have their counterparts 
in classical canonical variables, and thus provide a direct link between microscopic nuclear models and 
phenomenological collective models. A lot of attention has been paid' 14 ] especially after the success of 
the interacting boson model for nuclei J 15 l 

The purpose of this letter is to report an exact boson mapping of the reduced BCS pairing Hamiltonian. 
When m single fermions occupy the ji, ■ ■, jm levels, respectively, the reduced BCS pairing Hamiltonian 
for deformed nuclear system is given by 

m 

H B cs/G = $> t / G )( c It c ^ + cl L c jA ) + £ %ki - 'StSj, (1) 

T=l j i,j 

where Cj a (cj a ) are fermion creation (annihilation) operators, Sj~ = c jfC\j| and SJ — Cj^Cj^ are pair 

creation (annihilation) operators, kj — (cLc 3 -f + c^,Cji)/2, tj are single-particle energies taken from 
any deformed mean-field theory, G > is the equal strength pairing parameter, r\j = 2ej/G, and the 
summation sign with prime indicates that the sum is restricted to levels other than those occupied by 
the single fermions. 

Because solutions of m ^ cases are basically similar to those of seniority zero case, in the following 
we only consider the case with m = 0. For fc-particle excitation, the wavefunction of (1) in this case with 
p levels considered can be written as! 9, 10 1 

|fc; = S+(eP)S+(eP) ■ ■ ■ S+(E^)\0), (2) 
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where |0) is the pairing vacuum state satisfying Sj\0) = for 1 < j < p 



j=i V] - tip 



with the corresponing eigen-energy Ejp — G^*j =1 E^ 



The pair energies Ej£' should satisfy k coupled Bethe ansatz or Richardson-Gaudin equations: 

p i 2 

1 = E - T^y + E ( 4 ) 

for jit = 1, 2, • • • , k. It is understood that the additional quantum number £ in (2)-(4) is introduced to 
label the £-th set of roots {E^} of the equations (4). 

To map the reduced BCS Hamiltonian (1) into a boson Hamiltonian, we first use the mapping that 
maps the fermion pair operators kj, into the corresponding real boson operators with 

kj i ► rij = b^bj, Sf i-> &], Sy b 3 V j, (5) 

in which the images satisfy the usual commutation relations of boson operators: = and 

[bi,bj] = 0. It is clear that this mapping is different from that based on group structure' 14 ' because the 
images of no longer satisfy the commutation relations of the original SU(2) algebras. Furthermore, 
the mapping is unitary and number-conserving. We then seek a Bose Hamiltonian constructed from 
those boson images which should keep the wavefunction (3) consistent after the mapping. We found the 
one-body term in (1) does keep the same form after the mapping, while the pairing interaction term can 
not be mapped into one-body form, but with an additional non-Hermitian two-body interaction term, 
which is quite natural because the fermion pairing interaction like hard-core boson hopping can not be 
replaced by usual boson hopping. The final image of (1) after the mapping (5) is of the following form: 

H Bosc /G = ]T (ife - IK - E h \ h i + E n Mbj • (6) 
To reveal the dynamics of the Hamiltonian (6), let us consider a more general form of (6) with 

H B u = ^(2 £j - t - U)nj - X> - njU^bj + U £ n], (7) 

3 i=£j 3 

where 2e 3 - — t — U in the first term can be regarded as contribution from external potential or on-site 
disorder, the second term describes boson hopping among all sites with site dependent hopping parameter 
t—rijU, and the third term is the on-site repulsion. Since the two-body interaction term usually contribute 
with the same order of magnitude as the one-body term, the on-site repulsion parameter may be set as 
U = Uo/k, where k is the total number of bosons, in which Uq and t is of the same order of magnitude. 
Hence, the reduced BCS pairing Hamiltonian is mapped into a Bosc-Hubbard model with site-dependent 
hopping parameter t — Uo(rij/k). Therefore, the more the bosons on the j-th level, the less the hopping 
strength of other bosons hopping onto j-th level if t > J7o- I n the Bose Hamiltonian (6) with t = U = G, 
however, the condition 1 > nj is no longer satisfied if rij ^ or 1, which means that the fermion pairing 
interaction looks extremely repulsive after the boson mapping (5). 

To prove that (6) is indeed the exact boson image of (1), one can simply verify that the wavefunctions 
of (7), at least a part of them, can indeed be written as the boson image of (2) with 

\k,0 = B + (E^)B + (E^) ■ ■ ■ S+(4« ) )|0) B , (8) 
where |0)b is the corresponding boson vacuum, and 
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B + m=E - u ( 9 ) 



£J 2ej/t - 



with the corresponding eigen-energies 



and the Bcthe ansatz equations 



fc 



£«> = *5>«> (io) 



i = y 1 + y _™ (id 

for // = 1, 2, • • • , fc. 

Though it is difficult to analyze the spectrum generated by (11) analytically, one can verify that the 
whole spectrum described by (10) obtained from solutions of (11) are real and complete with (fc + p — 
l)!/(fc!(p — 1)!) eigenvalues when Uk/t < 1. Simple examples for p = 3, k = 2, 3 arc shown in Table 1. 

Table 1. The pair energies or roots of (11) and eigen-energies given by (10) for p = 3 and k = 2 and 3 with 
t = 1, ei = 1.1, e 2 = 2.2, £3 = 3.3, and U = 1/3. The dimension is exactly equal to (fc+p — l)!/((p— l)!fc!). 
k Dimension Eigenvalues Roots 



1.768 


£1 = 


0.884 - 0.792*' 


E 2 = 0.884 + 0.792z 






4.475 


Ei = 


0.935 


E 2 = 3.540 






6.763 


Ei = 


0.735 


E 2 = 6.028 






7.622 


Ei = 


3.811 - 0.358i 


E 2 = 3.811 + 0.358* 






9.744 


Ei = 


3.770 


E 2 = 5.974 






12.428 


Ei = 


6.214- 0.258z 


E 2 = 6.214 + 0.258* 






3.930 


Ei 


= 1.030 


E 2 = 1.450 - 1.229i 


£3 


= 1.450+ 1.229*- 


6.097 


Ei 


= 3.343 


E 2 = 1.377- 0.589* 


E3 


= 1.377 + 0.589* 


8.237 


Ei 


= 5.985 


E 2 = 1.126 - 0.685* 


E3 


= 1.126 + 0.685* 


8.577 


Ei 


= 1.227 


E 2 = 3.675 - 0.425*' 


E 3 


= 3.675 + 0.425*- 


10.704 


Ei 


= 3.672 


E 2 = 1.104 


E3 


= 5.928 


12.108 


Ei 


= 3.896 


E 2 = 4.106 - 0.568* 


E3 


= 4.106 + 0.568* 


13.302 


Ei 


= 0.932 


E 2 = 6.185 - 0.276i 


E 3 


= 6.185 + 0.276*- 


13.718 


Ei 


= 5.822 


E 2 = 3.948 - 0.302i 


E3 


= 3.948 + 0.302* 


16.142 


Ei 


= 3.848 


E 2 = 6.147- 0.304*- 


E3 


= 6.147 + 0.304* 


19.184 


Ei 


= 6.286 


E 2 = 6.449 - 0.357* 


E3 


= 6.449 - 0.357* 



When U = t = G, with which (7) is reduced to (6) corresponding to the image of the reduced BCS 
pairing Hamiltonian (1). In this case, though it is not guaranteed that all eigenvalues of (7) are real, 
especially for large k cases, a part of them satisfying (11) consisting of pl/((p — k)\k\) eigenvalues, which 
are the same as those given by (4), correspond exactly to those of the reduced BCS pairing Hamiltonian 
(1). There are [(fc + p — l)!/((p — l)!fc!) — pl/((p — k)\k\)] more other eigenvalues which are not provided 
by (10)-(11). Hence, eigenvectors other than those corresponding to the eigenvalues given by (10) can 
not be written in the Bethc ansatz form (8). 

Similar to the Dyson mapping! 14 ] and, for example, the iterative boson expansion approach, I 16 ! the 
resultant Bose Hamiltonian (6) is non-Hcrmitian. In addition, there will be spurious states involved in 
the boson space. Therefore, the Bose Hamiltonian should be projected onto the physical subspace. Let 
P be the projection operator, which can be expressed as 
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£ = 5>,0&*l, (12) 

where |fc, £) are normalized eigenvectors given by (8), and the sum runs over all possible values according 
to the number of solutions of (11) with U = t = G. Thus, it is clear that the projection operator P 
annihilates unphysical subspaceJ 14 ! It follows that the projected Bosc Hamiltonian with 

#Bosc = PHboscP (13) 

is digonalizable under the physical subspace spanned by {|fc, £)} with results shown by (8)-(ll). Hence, 
we obtain the Bose Hamiltonian i/eosc exactly equivalent to the reduced BCS Hamiltonian (1) in the 
physical boson subspace. 

In summary, an exact boson mapping of the reduced BCS pairing Hamiltonian is obtained under 
the guidance of the Richardson-Gaudin exact solutions of the reduced BCS pairing model. Though 
non-Hermitian, all solutions of the resultant Bose-Hubbard Hamiltonian are real and provided by the 
Richardson-Gaudin type wavefunctions and the corresponding Bcthe ansatz equations when Uk/t < 1. 
The physical Bose image of the reduced BCS pairing Hamiltonian is obtained by the projection method 
which is exactly equivalent to the reduced BCS Hamiltonian (1) in the physical boson subspace. Because 
what we have studied is based on deformed shell model like basis, the boson operators {bj,bj} do not 
conserve angular momentum which must be restored by angular momentum projection as mentioned in 
[17]. In the Nilsson mean-filed, for example, the boson operator b\ can be rewritten in terms of spherical 
pairs as 

b\ = 44 + 44 + x\gl + ■■■, (14) 

where x l L are transformation coefficients between the i-th Nilsson level and the spherical basis, and Sg , 
rfj, etc. are boson operators with L = 0, 2, • • •, and Ml = 0. After the angular momentum projection, 
one can better understand the intimate links between mean-field plus pairing models and the interacting 
boson model. Further study about this problem will be carried out in the near future. 

Support from the U.S. National Science Foundation (0500291), the Southeastern Universities Research 
Association, the Natural Science Foundation of China (10575047), and the LSU-LNNU joint research 
program (9961) is acknowledged. 

References 

[1] J. Bardeen, L. N. Cooper, and J. R. Schricffcr, Phys. Rev. 108, 1175 (1957). 
[2] M. Randeria, J. M. Duan, L.Y. Shieh, Phys. Rev. Lett. 62, 981 (1989). 
[3] D. J. Dean, Rev. Mod. Phys. 75, 607 (2003). 

[4] D. C. Ralph, C. T. Black, and M. Tinkham, Phys. Rev. Lett. 78, 4087 (1996). 

[5] L. Amico, A. Di Lorenzo, and A. Ostcrloh, Phys. Rev. Lett. 86, 5759 (2001). 

[6] D. W. Cooper, J. S. Batchlder, M. A. Taubcnblatt, J. Coll. Int. Sci. 144, 201 (1991). 

[7] P. Ring and P. Schuck, The Nuclear Many-Body Problem (Springer- Verlag, Berlin, 1980). 

[8] For example, J. Dukelsky, G. G. Dussel, J. G. Hirsch and P. Schuck, Nucl. Phys. A 714, 63 (2003), 
and references therein. 

[9] R. W. Richardson, Phys Lett. 3, 277 (1963). 



4 



[10] R. W. Richardson, J. Math. Phys. 6, 1034 (1965). 
[11] M. Gaudin, J. Physique 37, 1087(1976). 

[12] J. Dukelsky, C. Esebbag, and S. Pittcl, Phys. Rev. Lett. 88, 062501 (2002). 

[13] S. Rombouts, D. Van Neck, and J. Dukelsky Phys. Rev. C 69, 061303(R) (2004). 

[14] A. Klein, E. R. Marshalek, Rev. Mod. Phys. 63, 375(1991), and references therein. 

[15] F. Iachello and A. Arima, The Interacting Boson Model (Cambridge University Press, Cambridge, 
1987). 

[16] J. Dukelsky and S. Pittel, Phys. Rev. C 39, 2001 (1989). 

[17] F. Iachello and I. Talmi, Rev. Mod. Phys. 59, 339 (1987), and references therein. 



5 



